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Abstract 

We discuss local and global existence and uniqueness for the price formation free boundary model with ho- 
mogeneous Neumann boundary conditions introduced by Lasry & Lions in 2007. The results are based on a 
transformation of the problem to the heat equation with nonstandard boundary conditions. The free boundary 
becomes the zero level set of the solution of the heat equation. The transformation allows us to construct an 
explicit solution and discuss the behavior of the free boundary. Global existence can be verified under certain 
conditions on the free boundary and examples of non-existence are given. 

Resume 

Nous discutons l'existence locale et globale, ainsi que l'unicite des solutions pour le modele de formation des prix 
a frontiere libre avec des conditions aux bords de Neumann homogenes introduit par Lasry & Lions en 2007. Nos 
resultats sont bases sur une transformation de ce probleme en une equation de la chaleur avec des conditions 
aux bords non standard. La frontiere libre devient la ligne de niveau zero de la solution de l'equation de la 
chaleur. Cette transformation nous permet de construire une solution explicite et de discuter le comportement 
de la frontiere libre. L'existence globale peut etre verifiee sous certaines conditions sur la frontiere libre, et nous 
donnons des exemples de non-existence. 



1. Introduction 

In this paper we discuss the mathematical analysis of a price formation mean field model, as stated in 
[6] . The equation models the trading of an economic good between a large group of buyers and a large 
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group of vendors. It consists of a non-linear parabolic free boundary evolution equation which describes 
the evolution of the densities of buyers and vendors and also determines the price of the good. Here we 
consider the problem on the real interval {—L,L), where L denotes the maximum and — L the minimum 
price. The model is given by the free boundary value problem 

ft - f m = \(t)(5(x - p(t) + a) - S(x - p{t) - a)), x G (-L, L),t> 0, (la) 

A(t) = -/*&>(*),*). /(?(*),*)= 0, (lb) 

f(x, 0) = //, p(0) = po, for some po in (-L + a,L — a), (lc) 

f x (±L,t) = 0, t > 0, (Id) 

with < a < L and compatibility conditions at time t = 0: fi(po) = and fi(x) > for x < 
Po and fi(x) < for a; > po. Throughout this paper we assume that // is in L 2 (—L,L) and denote 
the L 2 norm on (—L,L) by ||-||. System (1) has been studied in a number of papers, see [3,7,2]. Here we 
present global (non-) existence results of a solution of (1) on the bounded interval (— L, L). 



2. Analysis of the Neumann problem - Transformation to the Heat Equation 

In the sequel we denote the positive and negative part of a function / defined for x G (—L,L) by 
/+ := max(/, 0), /~ := max(— /, 0). We extend f + ,f~ by the value outside the interval (— L, L). 
Proposition 2.1 Let f be a solution of the modified price formation problem consisting of equation (la) 
with conditions (lb), initial datum (lc) and modified boundary conditions (be): 

U-L,t)={° P (t)>-L + a /GM)= (0 pit)<L-a 

\f x (-L + a,t) p(t)<-L + a. xy \f x (L-a,t) p(t) > L — a. 

Then the function 



f + (x + na, t), x < pit) 
F(x,t) = { ^&=o_ V ^ W (3) 

■> / (x-na,t),x >p(t), 

* — ^n— 

is a solution of the following B VP for the heat equation 

F t = F xx , for all x G (-L, L), t > 0, (4a) 
F x (±L,t) = F x (±LTa,t), t>0 (4b) 
F(x,t = 0) =F r (x), for allx &{-L,L), (4c) 

where Fj is constructed from fj according to (3). Conversely, if F is a solution of (4) then f{x,t) — 
F(x,t) - F+(x + a,t) + F~(x - a,t) satisfies (la), (lb), (lc), (2). 

The construction of the function F is motivated by the fact that / + ( f~) has jump-discontinuities at 
x = pit) — a (x = p(t) + a) and x = p(t) of equal magnitude but opposite signs. The summation procedure 
then moves the jumps in the derivatives of / outwards. Note that the sum in (3) actually consists only 
of finitely many terms and that the free boundary p = p{t) of the price formation problem becomes the 
zero level set of the heat solution and vice versa. In principle the free boundary can leave and reenter the 
interval (— L + a, L — a) without impeding the existence of a global solution of (la), (lb), (lc), (2). 
We shall now construct an 'explicit' solution of (4) by separation of variables. We set F (x, t) = ip(x)ip(t) 
and find <p" /<p = ip/ifr = —z 1 . The boundary conditions give the equation G(z) := cos(zL)— cos(z(i— a)) = 
with the eigenfunctions (p(x) — Asin(zx) and H(z) := sin(zL) — sin(z(L — a)) = with eigen- 
functions (p = Bcos(zx). We easily compute G(z) = iff z = (27r/)/a, (2irl)/(2L — a) and H(z) — 
iff z = (2nl)/a,(Tr(2l — 1)/(2L — a) for / G Z. Next we conclude that H = H(z) is a sine-type 
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function of type L with simple and separated zeros. Therefore {exp(z^x), exp{i 2l~^r x )} * s a R-iesz 
basis in L 2 (—L,L), sec [1]. Also G = G{z) is a sine-type function of type L, its zero are separated 
except z = 0, which is a zero of order 2. Hence {exp(i^x), exp{i 2 2 ^ x), x}; g z is a Riesz basis in 
L 2 {~L,L), see Theorem D in [5]. By separately considering even and odd parts we conclude that 
{cos(^x), sin(^x), sin( 2 ^ a x), cos( ^^~^ x), 1, x}, I = 1,2,... is a Riesz basis of eigenfunctions of 
the heat equation (4). The solution of (4) is 

oo 

F{x, t) = J2[{Ai sm(uj lt ix) + Bi cos(wi j; x))e^i.<* + d sin^.j^e - ^.'* + D t cos^s^)^- 1 *] 

l=l ( 5 ) 
+A x + B , 

with 0J h i = (2nl)/a, Wa,i = {2ttI)/{2L - a), u 3 ,i = {{21 - l)n)/(2L - a). Note that A h B U C U D u A , B 
can be determined uniquely for every initial datum Fj E L 2 (—L, L). 

Theorem 2.1 (Global Existence) The BVP (la), (lb), (lc), (2) has a unique global solution f = 
f{x, t) for t > 0. Furthermore the free boundary p = p{t) is a smooth graph p{t) <G (-L, L) for all t > 0. 
The proof follows from the construction described above. The min-max principle implies that p = p{t) 
is a graph for t > and the Hopf principle implies that p is smooth. Finite-time oscillations of p(t) are 
excluded by the x-analyticity of solutions of the heat equations (see [2]). Since < J a F{x,t)dx and 

> J L _ F{x, t)dx are conserved in time (by the equation and the boundary conditions) we conclude 
that pttfe (-L, L) for all t > 0. 

Clearly, a solution of (la), (lb), (lc), (2) is a solution of (1) on a time interval [0, T] iSp(t) £ (— L+a, L—a) 
for t € [0, T}. Then the total mass of buyers Mb = J*^^ f(x, t)dx and vendors My = — f^. f(x, t)dx are 
time conserved quantities. 

Theorem 2.2 The BVP (1) has a global solution conserving the total mass of buyers and vendors iff the 
zero level set p of the solution of (4) satisfies p{t) 6 (-L + a, L — a) for all t > 0. Then the free boundary 
p(t) converges to p^ e (— L + a,L — a). 
The proof of Theorem 2.2 is based on the following lemmas. 

Lemma 2.3 The solution F converges exponentially fast to Fqo = A^x + B^ in L 2 {~L, L). 

Proof: From the Riesz base property we deduce that there exist Ci,C2 6 R + such that ci||Ff|| 2 < 

YZMi + B ? + C i + D l) + A 2 a + B 2 < ^IjF/ll 2 . For F = F — {A x + B ) we deduce 

dimi 2 < TZM" + Bf + °? + A 2 )]e- 2 ^ < e~ 2 ^ Y,7j A " + B ? + °l + D ") ± ^e-^UFzf, 

with the sequence ji = mm{{An 2 l 2 )/a 2 , (4tt 2 Z 2 )/(2L - a) 2 , {{21 - l) 2 7r 2 )/(2L - a) 2 ), I = 1,2,.... □ 
Lemma 2.4 The solution (5) satisfies \F x {x,t)\ < sup |(-FY)x| an d \F{x,t)\ < c, Vx e {—L,L), t > 0. 

x£(-L,L) 

Proof: The function V = F x satisfies the heat equation with V(—L, t) = V{—L+a, t), V{L, t) = V{L — a, t) 
and initial condition V{x,t = 0) = {Fj) x {x). If V assumes its maximum on the cylinder [— L, L] x [0, T) 
at either boundary x = ±L, then it must also assume a maximum (with the same value) in the interior 
(due to be on V). This contradicts the maximum principle, thus V must assume its maximum at t = 0. 
The same arguments hold for the minimum. 

Since F{x, t) = F(-L, t)+ f* L F x {y, t)dy we deduce aF{—L, t) = - f~^ +a F{x, t)dx+ f~^ +a f* L F x {y, t)dyd 
I\ + 12. We know that F x is bounded, therefore < K. In addition 4r J_^ +a F{x, t)dx = and there- 
fore J_L +a F{x,t)dx = f_L +a Fi{x)dx. Thus F{—L,t) as well as F = F{x,t) are bounded uniformly on 
{-L,L)xR+. □ 

Proof of Theorem 2.2: We know that F converges exponentially fast to F^ = A^x + B^, and that there 
exists a smooth graph p = p{t) such that F{p{t),t) = 0. Now we assume that p{t) 6 (— L + a, L — a) for 



t > 0, choose any sequence t n — > oo and conclude that there is a subsequence t nk such that p(t nk ) — > Poo £ 
[— L + a, L — a]. Let y> be a test function in T>(p oc ,L). If fc is sufficiently large we conclude f(x, t nk ) < 
for x £ supp<y9. Therefore f^ L f(x,t nk )tp(x)dx < 0. Since f(-,t) — > there is a subsequence i nfc( such 
that f(x,t nk ) converges to foo{x) pointwise a.e. in (—L,L). The function \f(x,t nk )\ < K on [— L,L\ 

for all fc. Then we deduce from Lebesgue's' dominated convergence theorem that J_ L f(x,t nk )<p(x)dx — > 

J^ L fooV{x) < 0. Since /oo = A 00 x + B 00 + (A 00 (x-a) + B 00 )~ -{A !yo (x + a) + B oc ) + we conclude < 
for X > poo and /oo > for x < Poo- From Lebesgue's' dominated convergence theorem, we conclude 
- fpu ) f(x, t)dx ->■ - foo{x)dx > 0. Analogously /coda > 0. 

Next we show that Foo has a unique zero in (-L, L). If i 7 ^ = AqoZ + i?oo > on (-L, L) then i 7 ]^ — 
0, i 7 ^ = -Foo and therefore /oo is constant in (-L, L), which is a contradiction. The same argument holds 
for Foe < 0. Therefore the function /oo is given by foo(x) = ±a for x £ (— £,Poo ~ a) and x £ (poo + ^) 
respectively and /oo = — a/a(x — p^) for x £ [poo — &,Poo + «], with a £ R + . We conclude that Poo is 
unique and that p(t) — > j>oo as t — > oo, since every sequence has a subsequence which converges to the 
same limit. □ 

Theorem 2.5 Let fj be such that Mb/Mv 4- [o/(4L — 3a), (4L — 3a)/o], where Mb, My denotes the 
initial mass of buyers and vendors. Then (1) does not have a global-in-time solution, which conserves 
both buyers and vendors masses. 

Proof: The result follows since a, Poo can not be adjusted such that poo £ [-L + a, L — a], where Mb = 
I-l foodx, M v = - foodx. □ 

Note that the choice of L, which corresponds to the maximally attainable price, is more or less arbitrary 
but a bad (too small) choice of L might impede global existence. In this case the model clearly looses 
its 'practical' significance. We remark that global existence results for (1) (with a free boundary which 
remains in {—L + a, L — a)) for initial data which are small perturbations of stationary solutions are 
straightforward, without using the analytical machinery of [4] . 
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